Macroscopic Einstein-Podolsky-Rosen pairs in superconducting circuits 



L.F. Wei, 1 ' 2 Yu-xi Liu, 1 Markus J. Storcz, 1 ' 3 and Franco Nori 14 

' Frontier Research System, The Institute of Physical and Chemical Research (RIKEN), Wako-shi, Saitama, 351-0198, Japan 
2 Institute of Quantum Optics and Quantum Information, Department of Physics, 
Shanghai Jiaotong University, Shanghai 200030, P.R. China 
3 Physics Department, ASC, and CeNS, Ludwig-Maximilians-Universitdt, Theresienstrasse 37, 80333 Miinchen, Germany 
4 Center for Theoretical Physics, Physics Department, CSCS, 
The University of Michigan, Ann Arbor, Michigan 48109-1040, USA 
(Dated: February 1, 2008) 

We propose an efficient approach to prepare Einstein-Podolsky-Rosen (EPR) pairs in currently existing 
Josephson nanocircuits with capacitive couplings. In these fixed coupling circuits, two-qubit logic gates could 
be easily implemented while, strictly speaking, single-qubit gates cannot be easily realized. For a known two- 
qubit state, conditional single-qubit operation could still be designed to evolve only the selected qubit and keep 
the other qubit unchanged; the rotation of the selected qubit depends on the state of the other one. These con- 
ditional single-qubit operations allow to deterministically generate the well-known Einstein-Podolsky-Rosen 
pairs, represented by EPR-Bell (or Bell) states. Quantum-state tomography is further proposed to experimen- 
tally confirm the generation of these states. The decays of the prepared EPR pairs are analyzed using numerical 
simulations. Possible application of the generated EPR pairs to test Bell's Inequality is also discussed. 

PACS number(s): 03.67.Mn, 03.65.Wj, 85.25.Dq. 



I. INTRODUCTION 

Quantum mechanics (QM) is a very successful theory. It 
has solved many physical mysteries in both macroscopic su- 
perconductivity and microscopic particles. Still, laboratory 
studies of its conceptual foundation and interpretation con- 
tinue to attract much attention. One of the most important 
examples is the well-known Einstein-Podolsky-Rosen (EPR) 
"paradox", concerning the completeness of QM. Based on a 
gedanken experiment, Einstein, Podolsky and Rosen (EPR) 
claimed |1] that QM is incomplete and that so-called "hid- 
den variables" should exist. This is because a two-particle 
quantum system might be prepared in a correlated (i.e., entan- 
gled) state, even though the two particles are spatially sepa- 
rated by a large distance and without any direct interaction. A 
measurement performed on one of the particles immediately 
changes the state (and thus the possible physical outcome) of 
the other particle. This "paradox" leads to much subsequent, 
and still on-going, researches. Bell proposed |2] an experi- 
mentally testable inequality to examine the existence of the 
hidden variables: if this inequality is violated, then there are 
no so-called local "hidden variables", and thus quantum me- 
chanical predication of existing quantum non-local correla- 
tions (i.e., entanglement) is sustained. 

During the past decades, a number of interesting experi- 
ments 1 3] using entangled photon pairs have been proposed 
and carried out to investigate the two-particle non-local corre- 
lations. These experiments showed that Bell's inequality (BI) 
could be strongly violated, and agreed with quantum mechan- 
ical predictions. Yet, one of the essential loopholes in these 
optical experiments is that the required EPR pairs were proba- 
bilistically generated in a small subset of all photons created in 
certain spontaneous processes. Thus, it is necessary to study 
two-particle entanglement in different, e.g., massive or macro- 
scopic systems, instead of fast-escaping photons. Expectably, 
the EPR pairs between these massive "particles" can be de- 



terministically prepared. Theoretical proposals include those 
with e.g., neutral Kaons |4], Rydberg atoms |5I], ballistic elec- 
trons in semiconductors |6], and trapped ions |7]. Experimen- 
tally, two Rydberg atoms had been first entangled to form EPR 
pair in a high Q cavity by the exchange of a single photon 1 8]. 
Later, by exchanging the quanta of the common vibrational 
mode, EPR correlations with ultralong lifetime (e.g., up to 5 
microsecond) had been generated between a pair of trapped 
cold ions 1 9 ] . Consequently, violations of BI have been exper- 
imentally verified with the EPR correlations between either 
the two ions iToll . or an atom and a photon fllll . 

Recent developments of quantum manipulation in coupled 
Josephson systems fL?! Fllll allow to experimentally investi- 
gate the quantum correlations between two macroscopic de- 
grees of freedom in a superconducting nano-electronic de- 
vice 0. Proposals have been made for producing quan- 
tum entanglement between two superconducting qubits, e.g., 
indirectly coupled by sequentially interacting with a current- 
biased information bus 0[kJ, coupled inductively fl7lfl8il . 
and coupled via either a cavity mode 1 19[ or a large Joseph- 
son junction 1 20]. By introducing an effective dynamical de- 
coupled approach, we have shown |21] that the BI could also 
be tested with superconducting qubits, even if the interaction 
between them is fixed. The robustness of the scheme pro- 
posed in Ref. l2~fll is better suited for weak interbit couplings, 
e.g., when the ratio of the interbit-coupling energy E m and the 
Josephson energy Ej of the qubit is small. In this paper, for 
an arbitrary interbit coupling strength, we discuss how to pre- 
pare the EPR correlations, i.e., deterministically generate and 
tomographically measure the well-known EPR-Bell (or Bell) 
states: 

|V±) = -^(|oo) ± |ii)), \<t>±) = -L(|oi) ± |io», (i) 

in a capacitively coupled Josephson circuit. Its possible appli- 
cation to directly test the EPR paradox is also discussed. 
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FIG. 1: (Color online) Two capacitively-coupled SQUID-based 
charge qubits. The quantum states of two Cooper-pair boxes (i.e., 
qubits) are manipulated by controlling the applied gate voltages 
V\ , V2 and external magnetic fluxes $1 , $2 (threading the SQUID 
loops). Pi and P2 (dashed line parts) read out the final qubit states. 



The outline of the paper is as follows. In Sec. II, a few ele- 
mentary quantum operations are proposed to deterministically 
manipulate two charge qubits coupled capacitively. Some of 
them only evolve a selected qubit and leave the remaining 
one unaffected. These operations are not strictly single-qubit 
gates (just conditional single-qubit operations), as the rotation 
of the selected qubit depends on the state of the other qubit. 
By making use of these operations, in Sec. Ill, we propose a 
two-step approach to deterministically generate the EPR pairs 
from the circuit's ground state \tp(0)) = |00). Further, we 
discuss how to experimentally confirm the generation of EPR 
pairs by tomographic measurements. In Sec. IV, considering 
the existence of typical voltage-noises and l//-noise, we nu- 
merically analyze the decays of the prepared EPR correlations 
within the Bloch-Redfield formalism 12211 . In Sec. V, we dis- 
cuss the possibility of testing BI with the generated EPR pairs. 
Conclusions and discussions are given in Sec. VI. 



II. MANIPULATIONS OF TWO CAPACITIVELY 
COUPLED JOSEPSHON CHARGE QUBITS 

We consider the two-qubit nano-circuit sketched in Fig. 1, 
which is similar to that in recent experiment JQEH- Two 
superconducting quantum interference device (SQUID) loops 
with controllable Josephson energies produce two Cooper- 
pair boxes, fabricated a small distance apart JT2II23I1 ) and cou- 
pled via the capacitance C m . The Hamiltonian of the circuit 
reads 



H 



E 

i=i,2 



E C] (hj 



"an 



-E { p cos 9j 



"9i> 



(2) 

in the charge basis. Here, the excess Cooper-pair number 
operator hj and phase operator dj in the jth box are conju- 
gate: [§j,hk] = iSjk- E C] = 4e 2 C Sfc /C s , j ^ k = 1,2 
and Ej — 2e,j j cos(7r$ : ,/$o) are the charging and Joseph- 
son energies of the jth box. E m — 4e 2 C m /Cs is the cou- 
pling energy between the boxes. Above, ej, and Ce, are 
the Josephson energy of the single-junction and the sum of 



all capacitances connected to the jth box, respectively. Also, 
Cs = C* Sl C* S2 - C 2 m and n 9] = C gj Vj/(2e). e is the elec- 
tron charge and $0 the A ux quantum. The circuit works in the 
charge regime with fc^T <C e,j j <C Ec, -C A, wherein 
quasi-particle tunnelling and excitation are effectively sup- 
pressed and the number rij (with nj = 0, 1, 2, ...) of Cooper- 
pairs in the jth boxe is a good quantum number. Here, 
fee, T, A, and 2ej. are the Boltzmann constant, temperature, 
superconducting gap, and maximal Josephson energies of the 
jth Cooper-pair box, respectively. 

Following Refs. II 2112311 . the dynamics of the system near 



the co-resonance point (where n gi 



' in 



= 1/2) can be ef- 



fectively restricted to the subspace S spanned by only the four 
lowest charge states: |00), 1 10) , |01) and |11), and thus the 
above Hamiltonian can be simplified to 



H = 



EM 



E12 o 



(3) 



with E 12 = E m /4, and E$ = E Cj (n gj - 1/2) + 
E m (n gk /2 — 1/4), j 7^ k = 1,2. The pesudospin opera- 
tors are defined as = |0j)(0j| - |lj){li| and = 
|0j)(lj| + |lj)(0j|. Here, the subindex j (or k) is introduced 
to label the state of the jth (or fcth) qubit. For example, |0j) 
refers to the logic state of the jth qubit is "0". For simplicity, 
the subindexes in a two-qubit state \mn) (with m, n = 0, 1) 
are omitted, and m (ri) usually (except when indicated other- 
wise) refers to the state |m) (|n)) of the first (second) qubit. 

Obviously, the interbit-coupling energy Eyi — E m /A is 
determined by the coupling capacitance C m and therefore is 

fixed by fabrication, i.e., not controllable. However, E^ and 

(7) 

Ej can be controlled by adjusting the applied gate-voltages 
Vj and fluxes $j, respectively. Although any evolution of this 
two-qubit system is solvable and can be expressed by a 4 x 4 
matrix in the subspace H, we prefer certain relatively simple 
quantum operations by properly setting the above controllable 
parameters to conveniently engineer arbitrary quantum states. 
These operations are summarized in the following three sub- 
sections. 



A. Operational delay 



First, we assume the circuit stays in the parameter settings 



such that E, 



c 



E 



U) 



0, until any operation is applied to 
it. Thus, during the operational delay r, the circuit evolves 
under the Hamiltonian H mt 
free time-evolution 



E12 (Tz^cTz , i.e., undergoes a 



U Q = 



{ 



V 



-ia 

















e lao 




Ei 2 

4^- (4) 



/ 



In this case, the Bell states in Eq. (1) will not evolve, once 
they have been generated. 
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B. Simultaneously evolving two qubits 



This implies that, if the fcth qubit is in the state |0fc) 



then the jth qubit undergoes a rotation A 



(A^P). During 



Due to the constant coupling, simultaneous operations on this operation the fcth qubit is unchanged and kept in its ini- 



two qubits are relatively easy. For example, if n 



1/2 (i.e., at co-resonance point) and E^p — E ( j' — Ej, then 



n(2) 



91 



' 92 



the circuit has the Hamiltonian H c , 



Ej{a x 1] +rt>)/2 + 



(2)s 



El2 which produces the following time-evolution 

operator 



E4o = 




(5) 



with 

a = cos(tQ/H) — iE\2 sin(iO/ft)/fi + exp(— itEi 2 /h), 
b = iEjsm(tn/h)/ft, Q= (Ej + Ef 2 )V 2 , 
c = cos(tn/H) — iEi2 sm(tfl/h)/n — exp(—itEi2/h). 

The subindex "co" refers to "co-resonance". Thus, we can 
simultaneously flip the two qubits, i.e., |00) ^ |11), and 
1 01) ^ 1 10} , by setting the duration as cos(if2/ft) = 
— cos(tEi2/h) = 1. Another specific two-qubit quantum op- 
eration 



fl-i l + i 

l + i 1-i 

1-i l+i 

V l + i 1-i 



(6) 



can also be implemented, if the duration is set as cos(ifi/ft) 
sm(tE 12 /h) = 1. 



C. Conditional rotations of a selected qubit 



Without the interaction free subspaces 1 24], a strict single- 
qubit gate cannot, in principle, be achieved in the system with 
strong fixed interbit-coupling. Recently, we have proposed 
an effective approach to approximately implement expected 
single-qubit logic operations 12111 . In what follows we show 
that conditional single-qubit operations, i.e., evolving only 
one selected qubit and leaving the other one unaffected, are 



(k) 



E 



(k) 







still possible. For example, one can set E { 
to only rotate the jth qubit. Indeed, the reduced Hamiltonian 



tial state. Obviously, if E i 
(thus cos aP 



U) 



2E\2 is satisfied beforehand 

0), and the duration is set as cos(tXj/h) = 

1, Xj = [(2E 12 ) 2 + (E ( f /2)] 1 / 2 , then the following two- 
qubit Deutsch gate Il25ll 



U®(8j) = I,- ® lOkXOjkl 



[I j cos 6j 



ia^' sin( 



\U)(U\, (8) 



with 9j — tEj/(2h), is obtained. Above, Ij is the unit op- 
erator relating to the jth qubit. The above operation implies 
that the target qubit (here it is the jth one) undergoes a quan- 
tum evolution, only if the control qubit (here, the fcth one) is 
in the logical state "1". If the duration is set to simultaneously 
satisfy the two conditions: sin (9j = 1 and cos(t\j/K) = 1, 
then the above two-qubit operation is equivalent to the well- 
known controlled-NOT (CNOT) gate, apart from a phase fac- 
tor. On the other hand, if Eq — —2E\ 2 is set beforehand, 
then the target qubit undergoes the same evolution only if the 
control qubit is in the logic state "0". The corresponding time- 
evolution operator reads 



U 



U), 



W = Jj®lU><i*l 

+ [/j-cos^+^sin^lOfeXOfcl. (9) 



Furthermore, if E 



(i) 



E. 



(2) 



E 



is set before- 



hand, then the above conditional operation (7) on the jth qubit 
(keeping the fcth one unchanged) reduces to 



U 



U) _ 



, 3 - B, ® |0 fe )(0 fe | +B] ® \l k )(l k \+^a 
with 

Q = cos(tjj/h) — i cosix,- sin(tjj/h), 



4,(10) 



i sin 



otj sa\(t^ij j ft) , cos a.j = Eyi 1 7j 



7j = ^{E 12 Y + {E { p/2f 
This operation can be further engineered to 



U U) -— 
Ui "V2 



(11) 



= E^'a^'/2 - E K j'a x v /2 + E 12 a^'a^' yields the if E { f = 2E 12 and sin. (7^/ ft) = 1 are further set. This is a 
following time-evolution 



U$ = A»> ® |0 fc )(0 fc | +A^> ® |l fc ><lfc|, 



50) 



(7) 



with 



.4 



0) 



^\o 1 )(o J \+^nw(h\+v± 1> ^>, 

lif = cos (t\ { l ] /ft) - i cosaf sin (tX { £> /ft) , 



,0) 



k 0) 



isirrn^ sin (U^/ftj , sina^ = E^'/(2X^'), 



>C0 



(i) 



'm>/2±E 12 } 2 



Hadamard-like operation on the jth qubit. 

Of course, the above operations, although only evolve the 
selected qubit and leave the other one unaffected, are no the 
strict single-qubit quantum gates (but just the especial two- 
qubit quantum operations). This is because the rotations of 
the selected qubit depend on the states of the other one. Note 
that, due to the presence of the constant interbit-coupling E12, 
(i) 

the value of E^ depends on both gate-voltages applied to the 

(2) 

two Cooper-pair boxes. For example, E c — requires that 
the two gate-voltages should be set to satisfy the condition: 

(n g2 -l/2)/(n gi -l/2) = -2E 12 /E C2 . 
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m. EPR-BELL STATES: THEIR GENERATIONS AND 
MEASUREMENTS 

Now, it will be shown how to deterministically generate 
EPR correlations between the above two capacitively coupled 
Josephson qubits. We will also propose how to experimentally 
confirm the expected EPR-Bell states. 

A. Deterministic preparations 

Naturally, we begin with the ground state of the circuit 
\ip(0)) = 1 00), which can be easily initialized by letting the 
circuit work far from the co-resonance point via a large volt- 
age bias. 

First, we prepare the superposition of two logical states of 
a selected qubit, e.g., the first one. This can be achieved by 
simply using a pulse of duration t\ to implement the above 
quantum operation (9), i.e., 

1-0(0)) = |00) |¥±) = -L(|00) ± i|10». (12) 

Here, the duration is set to satisfy the conditions 
cos(£iAi/?l) = landsin^i = ±l/v2. The plus sign cor- 
responds to the time durations for 6\ = 7r/4, and 3ir/4. The 
minus sign corresponds to 6i = 5tt/A, and 7tt/4. 

We next conditionally flip the second qubit, keeping the first 
one unchanged. The expected operations can be simply ex- 
pressed as either |00) — * 1 01), keeping 1 10) unchanged, or 
1 10) — > 1 1 1) , keeping |00) unchanged. The former (latter) op- 
eration requires to flip the second qubit if and only if the first 
qubit is in logic state "0" ("1")- These manipulations have 
been proposed above, and thus the desirable Bell states can be 
deterministically prepared by 

|* ± )^ 2 V±> = ^(|01)±|10)), (13) 

and 

l*±> + -VbA ± ) = ^(|00)±|ll)), (14) 

respectively. The duration t% of the second pulse is determined 
by the condition cos(A2t2/ft) = sin 62 = 1. 



B. Tomographic reconstructions 

The fidelity of the EPR correlations generated above can 
be experimentally measured by quantum-state tomography, a 
technique for reconstructing the density matrix of quantum 
state. For the complete characterization of an unknown two- 
qubit state with a 4 x 4 density matrix p — (pij,u) (with 
k,l = 0, 1), we need to determine 15 independent real 
parameters, due to tip = J2i,j=o.i PijM = and PvM = 
P%i ij ■ T ms can be achieved by a series of measurements on 



a sufficient number of identically prepared copies. The oper- 
ations presented above for the generation of EPR pairs could 
provide enough copies of any expected EPR pairs to be re- 
constructed. Experimentally, Bell states of pseudo-spins (e.g., 
in nuclear magnetic resonance systems 1 26], two-level trapped 
cold ions and the photon pairs l27ln have been tomograph- 
ically reconstructed by only using a series of single-qubit ma- 
nipulations. Recently, we have proposed a generic approach 
to tomographically measure solid-state qubits with switch- 
able interactions |28]. Due to the relatively strong interbit- 
coupling, which is always on in the circuits considered here, 
specific operations are required to realize the tomographic re- 
construction of the EPR pairs generated. 

The state of a charge qubit is often read out by capacitively 
coupling a single-electron transistor (SET) to the measured 
qubit 1 29]. When a projective measurement Pj = is 

performed on the state p, a dissipative current 1^ oc tr(pA ) 
flows through the jth SET coupled to the jth qubit. Such a 
projective measurement is equivalent to the measurement of 

Oz , as (Tz = (I — Pj)/2- For the present system one may 
perform three kinds of projective measurements: i) the P±- 
measurement (with projective operator P x ) acting only on the 
first qubit (independent of the state of the second qubit); ii) 
the P 2 -measurement (with projective operator P 2 ) operating 
only on the second qubit (independent of the state of the first 
qubit); and iii) the Pi 2 -measurement (with projective operator 
Pi <8> A) simultaneously acting on both Cooper-pair boxes. 

All diagonal elements of the density matrix p can be di- 
rectly determined by performing these three kinds of projec- 
tive measurements on the system. In fact, pi 1,11 can be deter- 
mined by the P12 -measurement as 

4 12) oc p n ,u = tr(pA <8 A). (15) 
Next, pio,io could be determined by Pi -measurement as 

IP oc pio.io + Pnai = tr(pA). (16) 
Also, we can determine poi.01 by the P2 -measurement as 

P? } oc p i,oi + pn,n = tr(pP 2 ). (17) 

The remaining element poo, 00 could be determined by the nor- 
malization condition trp = 1 . 

The 12 non-diagonal elements which are left, should be 
transformed to the diagonal positions of new density matrix 
p' = WpW\ by performing a proper quantum operation W 
on the original density matrix p. For example, after a quan- 
tum manipulation U^j \ see Eq. (11), evolving the system to 

p = U^P pt/j 1 ^, we can perform the P12 -measurement to 
obtain 

I { c 12) oc tr[pP®P 2 ] 

= 2^01,01 +Pn,ii - 2Re(poi,n)], (18) 

for determining Re(poi,n); an d perform the Pi -measurement 
to obtain 

7f cx tr[pA] = \[l +2Re(p o,io - Poi.u)], (19) 
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TABLE I: Tomographic characterization of an unknown two-qubit 
state p = (pij,ki) with i,j,k,l = 0, 1 in capacitively-coupled 
Josephson circuits. Each row of this table requires operating on an 
identically prepared initial state p. 



Operations 


Measurement 


Determining 


No 


Pl2 


/On, n 


No 


Pi 


(010,10 


No 


?2 


poi,oi 




Pl2 


Re(p i,n) 


&™ 


Pi 


Re(poo,io) 




Pl2 


Re(pio,n) 


up 


Pi 


Re(poo,oi) 


t r(l) f 7T\r"r(2) / 7T \ 


Pi 


Re(poo,n) 


Vfi^Vf (f) 


Pl2 


Re(poi,io) 




P 2 


Im(p o,io) 


t>| 15 (5) 


P 2 


Im(poi,ii) 


f>i 2) (z) 


P 2 


Im(/9 o,oi) 




P 2 


Im(pi ,ii) 


u co 


Pl2 


Im(p o,ii) 


u co 


P2 


Im(poi,io) 



for determining Re(poo,io)- All the remaining 10 off-diagonal 
elements of p can be similarly determined. 

Table I summarizes such a procedure for tomographic char- 
acterization of an unknown two-qubit state in this fixed- 
coupling two-qubit system. We need to first apply to p the 
quantum operations listed in the first column of Table I. After- 
wards, the projective measurements listed in the second col- 
umn of Table I must be made. In this way, all the matrix ele- 
ments of p can be determined. Of course, this is not a unique 
approach for determining all fifteen independent elements of 
the density matrix. In fact, the expected tomographic recon- 
struction could also be achieved by only using the Pi- and 
p2-measurements, and making the Pi2-measurement unnec- 
essary. 

With the density matrix p obtained by the above tomo- 
graphic measurements and comparing to the density matrix 
of ideal Bell states, i.e., 



P\1>± 



( 1 ±1 \ 



\±i o o i j 



P|0±> 



/ \ 
1 ±1 
±1 1 

V o o oo/ 



the fidelity of the EPR pairs generated above can be defined 



as F\ 



1>±) 



HPP\4>±)) andF | 



<t>±) 



tt (p/°i<*±))> res P ectivel y- 



So far, we have shown that EPR correlations could be pro- 
duced between two capacitively coupled Cooper-pair boxes. 
Further, these entangled states can be characterized by using 
tomographic techniques via a series of projective measure- 
ments. Below, we will numerically estimate the lifetimes of 
these states and discuss their possible application to test Bell's 
inequality. 



IV. DECAY OF EPR-BELL STATES DUE TO 
GATE- VOLTAGE NOISE 



The EPR pairs generated above are the eigenstates of the 

and thus are long-lived, at 



Hamiltonian H mt = E^ai 1 ^ ai 2 ^ 



least theoretically, in the idle circuit with Eq — 



Ef = 0. 



Under the influence of various disturbing perturbations, these 
pure quantum states will finally decay to the corresponding 
mixed states. In fact, experimental solid-state circuits are very 
sensitive to decoherence because of the coupling to the many 
degrees of freedom of the solid-state environment. However, 
coherent quantum manipulations on the generated EPR pairs 
are still possible if their decay times are sufficiently long. 



A. Model 

The typical dominating noise in Josephson circuits is 
caused either by linear fluctuations of the electromagnetic 
environment (e.g., circuitry and radiation noises) or by low- 
frequency noise due to fluctuations in various charge/current 
channels (e.g., the background charge and critical current fluc- 
tuations). Usually, the former one behaves as Ohmic dissipa- 
tion 13011 and the latter one produces a 1// spectrum |31], 
which is still not fully understood in solid-state circuits (see, 
e.g., jl3l). Here, we assume that the decay of the EPR pairs 
arises from linear environmental noises, i.e., we investigate 
the fluctuations of the gate voltages applied to the qubits. 
Moreover, the effect of background charges that cause dephas- 
ing are modeled by setting the zero-frequency part of the bath 
spectral function to a value given by the experimentally ob- 
tained 1 33] dephasing rates for the charge qubit system. This 
approach is valid for noise that can be approximated as lead- 
ing to an exponential decay. The effect of gate-voltage noise 
on a single charge qubit has been discussed in 13011 . We now 
study two such noises in a capacitively-coupled circuit. Each 
electromagnetic environment is treated as a quantum system 
with many degrees of freedom and modeled by a bath of har- 
monic oscillators. Furthermore, each of these oscillators is 
assumed to be weakly coupled to the Cooper-pair boxes. 

The Hamiltonian containing the fluctuations of the applied 
gate voltages can be generally written as 



with 



and 



H = H + H B + V, 



hb= [btfai + 5) 

3=1,2 Ulj V 7 



(20) 



V = (T { P{X X +(3X 2 )+ <rf ] {X 2 +7*1), (21) 

being the Hamiltonians of the two baths and their interactions 
with the two boxes. Here, 



E C 3 Cg 3 

4e 



Ylk&fiuj +gi 



(22) 
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with a LJj , being the Boson operators of the jth bath, and 
g UJj the coupling strength between the oscillator of frequency 
u>j and the non-dissipative system. Due to the mutual cou- 
pling of the two Cooper pair boxes, there will be crosstalk of 
the noise affecting each qubit. This is modelled in the spin- 
boson model with two bosonic baths represented above by the 
terms with the additional factors j3 and 7. The amount of 
this crosstalk is given by the network of capacitances or the 
corresponding energies only; namely, (3 — E m /2Ec 2 and 
7 = E m /2Ec 1 , and by inserting experimental values one 
finds that (3 w 7 » 1/10. 

The effects of these noises can be characterized by their 
power spectra. The spectral density of the voltage noise for 
Ohmic dissipation can be expressed as 

Jf(uj) = \g Uj \ 2 5(u>- Ljj) ~ r]huju? c l{u? c +ui 2 ). (23) 



Here, a Drude cutoff with cutoff frequency oj c = 10 4 GHz has 
been introduced, which is well above all relevant frequency 
scales of the system and given by the circuit properties 1 38 ] . 
The dimensionless constant 77 characterizes the strength of the 
environmental effects. Introducing the impedance, Z t (u>) = 
l/[iu>Ct + the spectral function for the fluctuations 

can be expressed via the environmental impedance Jf(ui) — 
ujRe(Zt(u))). Here, Z(uS) ~ Rv is the Ohmic resistor and Ct 
is the total capacitance connected to the Cooper-pair box. 

The well-established Bloch-Redfield formalism I22I l34ll 
provides a systematic way to obtain a generalized master 
equation for the reduced density matrix of the system, weakly 
influenced by dissipative environments. A subtle Markov ap- 
proximation is also made in this theory such that the result- 
ing master equation is local in time. In the regime of weak 
coupling to the bath and low temperatures, this theory is nu- 
merically equivalent to a full non-Markovian path-integral ap- 
proach I35I1 . For the present case, a set of master equations 
are obtained in the eigenbasis of the unperturbed Hamilto- 
nian 1 30] 



u nm ynrn 



kl 



Rnmki Pkl, 



(24) 



with the Redfield tensor elements R n mki given by 

r — A j- A V*r(~) _r(~) _r(+) 

*Himki — u £m 1 nrrk u nk L trrm L trunk 1 £mnfe' 



and the rates T^' given by the Golden Rule expressions 



(25) 



r (+) _ 
1 trunk " 11 



•.(-) 



dte-^ kt (V I/m (t)ViM0)), 



dt e 



'(ViMVViMt))- 



Here, Vij m (t) is the matrix element of the system-bath cou- 
pling term of the Hamiltonian in the interaction picture with 
respect to the bath, and the brackets denote thermal average. 



Note again that the strength of the dissipative effects is char- 
acterized by the dimensionless parameter r\. From experimen- 
tal measurements of the noise properties of the charge qubit 
system [36], it is found that the strength of the Ohmic noise is 
given by 



V : 



4e 2 i? 

f>7T 



1.8-10" 



(26) 



where R w 6 CI. Thus, current technology gives a noise floor 
of approximately 77 ~ 10~ 3 , which will be used for the nu- 
merical simulations. For visualization of the decay of the Bell 
states, we compute the concurrence 1 37], given by 



C 



c{0, 



g 4 }. 



(27) 



Here, the Qi, i — 1,2,3,4, are the eigenvalues of pp with 
p = (<Ty <g> <jy)p*(<7y ® ct 2 ,). The concurrence is a measure 
for entanglement and indicates non-locality. The maximally 
entangled Bell states (i.e., the ideal EPR correlations) yield a 
value of 1, whereas a fully separable state gives 0. 



B. Numerical results 

The results of the simulations are shown in Fig. [2] where 
the time evolution of the concurrence C shows the decays of all 
Bell states, for temperature set to an experimentally feasible 
value of 10 mK. The lifetimes of the operationally idle EPR 
pairs are of the order of several p s and thus sufficiently long 
(compared to the duration ~ 100 ps of the usual quantum 
manipulation). 

For the case where only the coupling term between the 
qubits is present and all single-qubit terms in the Hamilto- 
nian are suppressed, Fig. 2(a) shows that the Bell states 
decay exponentially fast to zero: C(t) ~ exp(— At), with 
A ~ 2.13 x 10 6 for \<f>±) and A ~ 3.18 x 10 6 for |^±). In 
this case, only pure dephasing contributes to overall decoher- 

ence rates, as H = H int = E^ 2 &z <Tz and [H, V] = 0, 
see Ref. 13811 . The magnitude of the dephasing part of deco- 
herence is essentially determined by the 1/f-noise. To model 
this, a peak in the spectral function at zero frequency can be 
introduced with a magnitude given by microscopic calcula- 
tions or experimental measurements of the magnitude of 1//- 
noise in these qubit structures. However, note that often the 
noise leads to non-exponential decay, which can neither be 
modeled by Bloch-Redfield theory nor be parametrized by a 
single rate. Here, we assume Markovian and Gaussian noise 
and set the zero frequency contribution, i.e., the dephasing 
due to the 1 / /-noise to an experimentally reported value of 
T v ss 10 7 Hz 1 33]. Note that the individual contributions 
from different noise sources sum up in the spectral function 
Jj:(oj) = Jf{oj) + J\/f(oS), which also holds at oj = 0. It 
is interesting to note that the decay time is independent of the 
inter-qubit coupling strength E12. In more detail, when the 
coupling energy E12 in the Hamiltonian is increased the de- 
cay does not change. The reason for this behavior is that the 
pure dephasing is only affected by the zero frequency part of 
the spectrum, which is obviously independent of the individ- 
ual frequency splittings, i.e., the characteristic energy scale 
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(a) 



ly >, ly > 




0.4 0.6 

t(ns) 



FIG. 2: Simulated time evolution of the concurrence C for a two- 
qubit system coupled to a noisy environment and initially prepared 
in the Bell states. Here, the temperature and the strength of noise are 
set to T = 10 mK and r\ = 10~ 3 , respectively, (a) captures the long- 
time decay of the concurrence for different entangled input states 
in the case of vanishing single-qubit terms, i.e., when only the inter- 
qubit coupling terms are present, (b) compares the decays of \ip-) for 



different interbit-couplings (E m = Ej, and O.lEj) without (E 



(i) 



E^p = 0), and with Josephson tunneling (E^ 
55 /xeV). 



E 



(2) 



Ej 



of the Hamiltonian. Also, one of the most important results, 
namely that the decay time of \<t>±) is longer than that of | ?/>+), 
is consistent with the analog experimental one in ion traps |9]. 
This is because \4>±) is the superposition of the two states with 
the same energy, while \ip±) corresponds to higher energy and 
is more sensitive to such perturbations. 

When the Josephson-tunneling terms exist, e.g., Ej — 

(2) 

Ej = Ej, we see from Fig. 2(b) that the decays of the 
generated EPR pairs are significantly faster than in the for- 
mer case without any tunneling. This is becasue the addi- 
tional Josephson tunneling provides additional decoherence 
channels since the Hamiltonian of the circuit now does not 
commute with the couplings to the baths. Moreover, also the 
overall energy scale in the Hamiltonian increases. In this case, 
the weaker interbit-coupling corresponds to the slower decay 
of the EPR pairs. 



V. TESTING BELL'S INEQUALITY 

A possible application of the deterministically generated 
EPR pairs is to test BI at the macroscopic level. Due to the 
existence of interbit constant-coupling, the required local op- 
erations of encoding classical information {8j} into the EPR 
pairs cannot be strictly implemented. In Ref. l2~fll we pro- 
posed an approach to overcome this difficulty by introduc- 
ing the effective single-qubit operations including corrections 
due to the constant-coupling. Instead, here we approximately 
perform the encoding procedure by sequentially applying the 



conditional single-qubit operations Uj , (j = 1,2) in Eq. 
(10). For the case of a\ — a 2 — a, the validity of the 
above quasi-local encodings could be described by the vari- 
ation of the degree of entanglement (i.e., concurrence) of the 
EPR pairs 



AC = 1 - y/l - [sin(2a)(l - cos(2</>i + 2tp 2 ))/2} 2 , (28) 

with ifj = 2jjt/h. Obviously, AC — corresponds to the 
ideal locality or maximal locality. After the above encoding, 
we simultaneously detect [13] the populations of qubits and 
check if they are in the same logic states: the excited one |1) 
or the ground state |0). 

Theoretically, the correlation of two local variables, ipx and 
if2, can be defined as the expectation value of the operator 



P T = |11)<11| 
and reads 



- 100)(00| - |10)(10| - |01)(01| = <g> 



•(2) 



E(<Pl,<P2) 



cos 2 a 



■ sin 2 a cos(</?i + ip 2 ) 



(29) 



Experimentally, all the above operational steps can be re- 
peated many times in a controllable way for various param- 
eter sets. As a consequence, the correlation function E can be 
measured by 



N aa , me ((px,ip 2 ) +Ndis(<pi,ip2)'' 



(30) 



for any pair of chosen classical variables ipi and (f2. Here, 
N same ((px,cp2) (N dit f(ipi,ip 2 )) are the number of events with 
two qubits found in the same (different) logic states. With 
these measured correlation functions, one can experimentally 
test the BI in the present superconducting systems. 

We consider the following typical set of angles: 
{<Pj,<p'j} — {— 7r/8, 37r/8} and the interbit couplings 
E m = 4E 12 = Ej,Ej/10, and £j/100, respec- 
tively. The corresponding variations AC of the concur- 
rence and the correlation E(<pi,ip 2 ), which yields the 
Clauser, Home, Shimony and Holt (CHSH) (3J function 
/ = \E^ U ip 2 ) + E{^i P2 ) + E{ip U ip'. 2 )-E{ip' 1 ^' 2 )\, 
are given in Table II. It is seen that the variations AC of 
the concurrence, after the above quasi-local operations 



U 



0) 



decrease with decreasing interbit coupling. For very 



weak coupling, e.g., E m /Ej — 0.1(or0.01), the applied 
conditional single-qubit operations can be regarded as local, 
away from 0.4%, (or 0.004%). Besides these tiny loop- 
holes of locality, Table II shows that the CHSH-type Bell's 
inequality |3J 



/(lv>V» < 2 



(31) 



is obviously violated. 



VI. DISCUSSION AND CONCLUSION 

Similar to other theoretical schemes (see, e.g., Ref. |!l||]) the 
realizability of the present proposal also faces certain techno- 
logical challenges, such as the rapid switching of the charge- 
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TABLE II: Variations of the concurrence, AC, correlations E, and 
CHSH-functions /, for certain typical parameters of the interbit cou- 
pling E m and the controllable classical variables ipi and if 2- 



E m 


(ipi, 1P2) 


AC 


E((pi, 1P2) 


/ 


Ej 


(-—71-/8 -tt/8) 


0.00699 


0.76569 


2.6627 


(— 7r/8,37r/8) 


0.00699 


0.76569 


(3tt/8, -tt/8) 


0.00699 


0.76569 


(37r/8,37r/8) 


0.26943 


-0.36569 


Ej/W 


(-tt/8, -tt/8) 


0.00238 


0.72434 


2.8264 


(-tt/8, 3tt/8) 


0.00011 


0.70784 


(3tt/8, -tt/8) 


0.00011 


0.70784 


(37r/8,37r/8) 


0.00363 


-0.70285 


Ej/WO 


(-tt/8, -tt/8) 


0.00001 


0.70711 


2.8284 


(-71-/8,371-/8) 


0.00001 


0.70711 


(3tt/8, -tt/8) 


0.00001 


0.70711 


(37r/8,37r/8) 


0.00004 


-0.70706 



and Josephson energies of the SQUID-based qubits and deco- 
herence due to the various environmental noises. Our numer- 
ical results, considering various typical fluctuations, showed 
that the lifetime of the generated EPR pairs adequately allows 
to perform the required operations for experimentally testing 
Bell's inequality. Indeed, for current experiments 1 12], the de- 
cay time of a fwo-qubit excited state is as long as ~ 0.6 ns, 
even for the very strong interbit coupling, e.g., E m ~ Ej. 
Longer decoherence times are possible for weaker interbit 
couplings. In addition, for testing this, the influence of the 
environmental noises and operational imperfections is not fa- 
tal, as the nonlocal correlation E(ipi, ipj) in Bell's inequality 
is statistical — its fluctuation could be effectively suppressed 
by the averages of many repeatable experiments. 



In summary, for the experimentally realized capacitively 
coupled Josephson nanocircuits, we found that several typi- 
cal two-qubit quantum operations (including simultaneously 
flipping the two qubits and only evolving a selected qubit in 
the case of leaving the other one unchanged) could be easily 
implemented by properly setting the controllable parameters 
of circuits, e.g., the applied gate voltages and external fluxes. 
As a consequence of this, macroscopic EPR correlated pairs 
could be deterministically generated from the ground state 
1 00) by two conditional single-qubit operations: prepare the 
superposition of the two logic states of a selected qubit, and 
then only flip one of the two qubits. To experimentally con- 
firm the proposed generation schemes, we also propose an ef- 
fective tomographic technique for determining all density ma- 
trix elements of the prepared states by a series of quantum pro- 
jective measurements. The deterministically generated EPR 
pairs provide an effective platform to test, at the macroscopic 
level, certain fundamental principles, e.g., the non-locality of 
quantum entanglement via violating the Bell's inequality. 

The approach proposed here can be easily modified to engi- 
neer quantum entanglement in other "fixed-interaction" solid- 
state systems, e.g., capacitively (inductively) coupled Joseph- 
son phase (flux) system and Ising (Heisenberg)-spin chains. 
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